Quantum calibration of measuring apparatuses 



^ ■ 

O ■ 

o : 

^ ■ 
<: 

00 , 



> 



oo : 
o . 

o • 

Oh' 



Giacomo Mauro D'ArianoQ Lorenzo Maccone0 and Paoloplacido Lo Prest^ 
QUIT - Quantum Information Theory Group, Dipartimento di Fisica 
"A. Volta" Universita di Pavia, via A. Bassi 6, 1-27100 Pavia, Italy. 
(Dated: February 1, 2008) 

By quantum calibration we name an experimental procedure apt to completely characterize an 
unknown measurement apparatus by comparing it with other calibrated apparatuses. Here we show 
how to achieve the calibration of an arbitrary measuring apparatus, by using it in conjunction with 
a "tomographer" in a correlation setup with an input bipartite system. The method is robust to 
imperfections of the tomographer, and works for practically any input state of the bipartite system. 
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The calibration of measuring apparatuses is at the ba- 
sis of any experiment. Theory and experiment are un- 
avoidably interwoven, and the calibration procedure of- 
ten needs a detailed knowledge of the inner working of the 
apparatus, especially at extreme precisions and sensitiv- 
ities, where a quantum mechanical description is needed. 
Here, the actual "observable" that is measured depends 
crucially on the microscopic details of the apparatus, and 
without knowing them the measurement lacks physical 
interpretation. 

In a quantum mechanical description, the calibration 
of a measuring apparatus corresponds to the knowl- 
edge of its POVM (positive operator-valued measure 0), 
which gives the probability p{n) of any measurement out- 
come n for arbitrary input state, via the Born rule 



p{n) Tr[pP„]. 



(1) 



In Eq. p is the density operator of the state on the 
Hilbert space Ti of the system, and the POVM is given 
by the set of operators {Pn} on H. To ensure that p{n) is 
a probability, the POVM must satisfy the positivity and 
normalization constraints P„ ^ 0, ^„ P„ = /. 

The concept of POVM generalizes the familiar von 
Neumann observable describing perfect measurements. 
Here the probability of obtaining outcome n is given by 
p{n) = K-^jon)!^, {|o„)} denoting a complete orthonor- 
mal basis for H, i. e. with POVM given by the one- 
dimensional projectors P„ — |o„)(o„|. The physical in- 
terpretation of the measurement is given via a quanti- 
zation rule that associates a self-adjoint operator O to a 
classical observable, |o„) being the eigenvector of O corre- 
sponding to its nth eigenvalue o„. However, this concept 
of observable does not cover many practical situations — 
e. g. phase-estimation|^ ISl, joint measurements of in- 
compatible observables |j, |S| , discrimination among non- 
orthogonal states 0, , informationally complete mea- 
surements 0, transmission of reference frames — and 
here the POVM description is needed. But then, in ab- 
sence of a direct physical interpretation of the measure- 
ment, we are faced with the problem of assessing the 
correct functioning of the measuring apparatus. 

Inferring the POVM of an apparatus through the the- 



oretical description of its functioning leads to quite in- 
volved derivations, based on different kinds of approx- 
imations. A paradigmatic case is that of the photo- 
counter [loj , where the number of photons claimed to be 
detected — usually very uncertain — is typically inferred 
from the cascading mechanism of the amplification pro- 
cess. The calibration is given essentially in terms of 
quantum efficiency and dark-current, and mostly satu- 
ration effects categorize detectors into the major classes 
of "linear" and "single-photon" . Even in a very simpli- 
fied model, a theoretical description accounting for the 
above features is very involved jll. ,12.] , and the resulting 
theoretical calibration is exceedingly indirect. 

The above scenario raises the following problem: is 
it possible to calibrate a measuring apparatus — i. e. to 
determine its POVM — with a purely experimental proce- 
dure, e. g. by comparing the apparatus with other (previ- 
ously calibrated) apparatuses? In this paper we propose 
a method to determine a POVM experimentally. The 
method uses the unknown apparatus jointly with a cal- 
ibrated "tomographer" on a suitably prepared bipartite 
system, as in Fig. ^ and the calibration results from the 
analysis of the correlations of outcomes. [A tomographer 
is an apparatus that measures an observable tunable in 
a complete set called quorum: more details on quantum 
tomography will be given in the following]. The basic 
scheme of the method stems on a previous method for the 
tomographic reconstruction of quantum operations 1^, 
and generalizes a popular calibration scheme de- 
signed to determine the quantum efficiency of a photo- 
detector. As it will be shown in the following, there is 
ample freedom in the choice of both the input bipartite 
state and the tomographer. The joint measurement must 
be repeated many times, analyzing the measurement out- 
comes with a proper tomographic algorithm [Hill, the 
POVM calibration being approached in the limit of in- 
finitely many outcomes. For finite set of data, the re- 
constructed POVM will be affected by statistical errors, 
which can be precisely estimated via the tomographic 
algorithm. The method works for generally infinite- 
dimensional Hilbert space (yielding a finite number of 
POVM elements, corresponding to the actually occurred 
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outcomes). 





FIG. 1: Experimental setup to determine the POVM of the 
unknown measurement apparatus A. The apparatus A is used 
jointly with a "tomographer" on a bipartite system prepared 
in a pre-determined state R. The tomographer measures an 
observable B^''^ from the quorum {B^''^}, yielding result m, 
whereas the unknown apparatus gives outcome n. The joint 
outcomes (n, m) are then processed using a tomographic al- 
gorithm, to finally obtain the POVM {P„} of A. 

The following simple example illustrates how the pro- 
cedure works. Suppose we know that the apparatus mea- 
sures an observable, but we don't know which one, and 
denote it by the orthonormal basis {|o„)}. We can use the 
maximally entangled input state — X)iLi N)N)/v^ 
the space Ti®T,T denoting the space of the quantum 
system impinging into the tomographer. The state can 
be equivalently written as |*) = ^ \oj)\o*) (|o*) 

denotes the vector with the complex conjugated coeffi- 
cients of \oj) with respect to the basis Then, the 
outcome n of the unknown measuring apparatus condi- 
tions the state pn = |o*)(o*| at the tomographer, and 
the POVM can be recovered using state reconstruction. 

We now present the general quantum calibration pro- 
cedure. Let's fix one observable at the tomographer, and 
denote it by {|6„)}. Upon denoting by {P„} the POVM 
of our measuring apparatus that we want to calibrate, 
the Born rule predicts that the outcome (n, m) of the 
joint measurement will occur with probability 



p(n,m) =Tr[(F„® |6„,)(6™|)i?]. 



(2) 



where R is the joint state of the two quantum systems, 
and we remind that the POVM of the joint measure- 
ment is given by the tensor product of the individual 
POVM's. Upon rewriting the joint probability in terms 
of the conditional probability p{m\n) via Bayes' rule, we 
conveniently introduce the state p„ at the tomographer 
conditioned by outcome n at the unknown measuring ap- 
paratus, namely 

p(n,m) =p{n)p{m\n) = p(n)Tr [p„|6„i)(6,„|] . (3) 

Upon evaluating the trace in Eq. (01 m two steps, i. e. 



and by equating Eqs. Q and Q for any possible vector 
\bm) (i- e. any possible observable), we have Pnp{n) = 
Tri[(P„ (g) IL)i?], namely 

The POVM element P„ can be recovered from the con- 
ditioned state Pn as follows 



Pn=p{n)n-\pn), 
by inverting the map 

n{X) = Tri[(X® ]l)i?]. 



(6) 



(7) 



p{n,m) = Tr |6™)(5,„|Tri [(P„ ® t) R] 



(4) 



X denoting an operator on Ti. The map TZ depends only 
on the input state i?, which then must be known. Hence 
we need a pre- calibration stage in which we previously 
determine the joint state R (this can be done via a joint 
quantum tomography with two equal tomographers on 
the input state R). Invertibility of the map TZ corre- 
sponds to a so-called faithful state |l3l |. Since invertible 
maps are a dense set, then almost any quantum state R 
is faithful. Of course, when approaching a state corre- 
sponding to a non invertible map, some information on 
the POVM {Pn] will be lost, corresponding to increas- 
ingly large statistical errors for some matrix elements of 
the operators P„ (inverting a linear map is clearly equiv- 
alent to inverting an operator: for the reader who prefers 
operators to map, an explicit connection between opera- 
tors and maps is given in Ref. 0|). 

Once the inverse map Hr^ has been calculated, we 
use quantum tomography in order to recover p„. Shortly 
quantum tomography is a method that allows us to es- 
timate the ensemble average of an arbitrary (complex) 
operator X by measuring a set of observables {P'^'^-'}, 
called quorum, which span the space of operators of the 
systern (for recent reviews on quantum tomography, see 
Refs. 0, 01). Typical examples of quorums are the 
three Pauli matrices Ux, Ty, for a qubit, or the set of 
quadratures = ^(0^6*"^ -I- ae~*"^) for a single mode of 
the radiation field with annihilation and creation opera- 
tors a and a\ respectively, g [0,7r) playing the role 
of the observable label within the quorum {X^f,} [X^ 
is measured by a homodyne detector at phase 4> rela- 
tive to the local oscillator 10). In short, the generic 
operator X is expanded as X ^ Xlfe Tr[XC('=)1']P('=), 
{C^*^^} denoting a dual set of the quorum {P^*^^} (the 
sum being replaced by an integral for continuous k), the 
quantity Tr[XC''^^^] being evaluated analytically. No- 
tice that one can remove the effects of noise at the to- 
mographer if the noise map M is invertible, by writing 
X = Y.k Tr[AA-i(X)C('=)t]AA(p(fe)). 

For the tomographic reconstruction we can either: a) 
average over the quorum, e. g. estimate {X) via the 
ensemble averages of the quorum observable as {X) = 
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X;feTr[XC('=)'f](S('=)) (the estimation of the density ma- 
trix element pij corresponding io X = b) we 
can use the maximum hkehhood approach |l7| . In this 
case, the estimated POVM elements P„ will maximize 
the probability Tr[(P„ (g) \bm^ ) {bm'' \) R] in the joint mea- 
surement on the pre-determined input state R of getting 
outcome n on the unknown measuring apparatus and m 
for the fcth observable i?'-'^-' of the quorum. Equivalently, 
one can maximize the logarithm of this quantity and con- 
sider simultaneously all the N joint measurement out- 
comes, corresponding to maximizing the likelihood func- 
tional 



N 



(8) 



under the constraints P„ ^ and J^n-^^ = I- Other 
prior knowledge about P„ can be easily incorporated 
by adding further constraints. Moreover, we can ac- 
count for a known source of noise Af at the tomographer, 
by replacing the projector \bml'){bm'i^\ in Eq. ((SJ with 

m\bti'){bl^i^\). 

Therefore, the procedure to calibrate an unknown mea- 
surement apparatus can be summarized in the following 
steps: i) [pre- calibration] using two tomographers, recon- 
struct the input joint state R. Check whether R is faith- 
ful; a) \joint measurements with the unknown apparatus] 
Replace one tomographer with the unknown detector, 
and collect N pairs of outcomes {ni,mi\, i — 1,...,N in 
a set of joint measurements with randomly selected ob- 
servable iJ^*^') in the quorum; Hi) [Data analysis] From 
the experimental data collect the probability p{n) of 
the outcome n at the unknown measurement apparatus, 
and then estimate the POVM {P„} using a tomographic 
strategy — either the averaging or the maximum likeli- 
hood. In the first case evaluate the density matrix pn 
of the state impinging in the unknown measuring appa- 
ratus, and then use Eq. (jBJ to recover the POVM. In the 
second case, evaluate the POVM directly by maximizing 
the likelihood functional £ in Eq. (jHJ on the given set of 
experimental data, with the state R obtained at step i). 

In Fig. 13 we present a simulated experiment of the 
quantum calibration of a photo-counter using homodyne 
tomography with the averaging strategy. The model 
of the calibrated detector is given in Fig. |31 Since 
the resulting POVM is diagonal in the photon-number 
basis, we limit the reconstruction to the diagonal ele- 
ments only. As input state R we use a twin beam state 
from parametric down-conversion of vacuum, of the form 
(X J2m^"^\^) ® I™)' where ^ is related to the amphfi- 
cation gain and |m) denotes the eigenstate of the pho- 
ton number. One can easily check that the twin beam is 
faithful for all ^ 7^ 0. As typical imperfection of the tomo- 
grapher, we consider non-unit quantum efhciency rjh for 
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FIG. 2: Calibration of the photodetector in Fig. |3 with 
rjp = 80% and u — 1, using a twin-beam with ^ = 0.88 (see 
text), and homodyne tomography with quantum efficiency 
rjh ~ 90%. The plots are the reconstruction of the diagonal 
matrix elements {n\Pk\n) of the detector POVM, from a set of 
5 X 10® compute-simulated data, using the averaging strategy. 
The reconstructed POVM is at the middle of the error-bars, 
whereas the theoretical values, for comparison, are given by 
the thick line. 
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FIG. 3: Model of the photodetector calibrated in the simu- 
lations of Figs. 121 and 2] Non-unit quantum efficiency rip and 
dark-current with averaged photon number v are equivalent 
to preceding an ideal detector by a beam-splitter of transmis- 
sivity rip mixing the input signal with a thermal mode with f 
average photons. 



the homodyne detector (the noise map can be inverted as 
long as 77?! > i 01). Since we reconstruct only the diag- 
onal part of the POVM, one can easily show that there 
is non need of knowing the homodyne phase 0, which, 
however, must be randomly distributed (the knowledge 
of (j) would allow to recover also the off-diagonal elements 
of the POVM). 

In Fig. ^ we present the same calibration, but using 
the maximum likelihood strategy. The convergence of 
the maximum-search algorithm is assured by the strict 
convexity of the likelihood functional £ over the space 
of diagonal POVM's (the convergence speed, however. 
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FIG. 4: The same as in Fig.|21 but using the maximum hkeli- 
hood method. Here only 5 x 10* simulated data are used. The 
error bars are obtained by standard bootstrapping techniques 
over a virtual repetition of 50 experiments. Notice how the 
result is statistically less noisy than that in Fig. |5| even for a 
10~^ smaller set of data. 



can be practically very slow). In the simulation we used 
a blend of sequential quadratic programming routines 
(to perform the constrained maximization) along with 
expectation-maximization techniques [l^ . By compar- 
ing Figs. |21and01we can see how the maximum likelihood 
estimation is more statistically efficient (i. e. fewer data 
are needed to achieve the same statistical error) than the 
averaging strategy and, in addition, the maximiza- 
tion of the likelihood recovers all the POVM elements 
simultaneously. On the other hand, compared to the av- 
eraging strategy, the maximum likelihood approach has 
the drawback of being biased, since one needs to put a 
cut-off to the Hilbert space dimension of the tomographic 
reconstruction and/or to the cardinality of the POVM. 
Both simulated experiments use realistic parameters and 
are feasible in the lab with current technology (see, for 
example, Refs. IH IH HHl^ ) , the major challenge of 
a real experiment being the matching of modes between 
photo-counter and homodyne detector, also ensuring that 
the detected modes are the same of the pre-calibration 
stage. 

We acknowledge financial support by INFM PRA- 
2002-CLON and MIUR for Cofinanziamento 2003. PL 
acknowledges partial support from ATESIT project IST- 
2000-29681. 



* Electronic address: |dajianoQunipv.rt| Also at Center for 
Photonic Communication and Computing, Department 
of Electrical and Computer Engineering, Northwestern 
University, Evanston, IL 60208 



[1 
[2: 

[3; 
[4: 
[s: 

[6: 

[7: 
[s: 

[9: 

[10 

[11 

[12 
[13 
[14 
[15 
[16 

[17 
[18 

[19 

[20 
[21 
[22; 
[23; 



maccone 


Sunipv.itl 


lopresti(§ 


unipv.iv 



C. W. Helstrom, Quantum Detection and Estimation 
Theory (Academic Press, New York, 1976). 

Physica Scripta T48 (1993) (special issue on Quantum 
Phase and Phase Dependent Measurements) . 

G. M. D'Ariano, C. Macchiavello and M. F. Sacchi, Phys 
Lett. A (248) 103-108 (1998) 

E. Arthurs and J. L. Kelly, Bell. Syst. Tech. J., 44 725- 
729 (1965). 

J. P. Gordon and W. H. Louisell, in Physics of Quan- 
tum Electronics, pp. 833-840, McGraw-HiU, (New York, 
1966). 

H. P. Yuen, R. S. Kennedy, and M. Lax, IEEE Trans. 
Inf. Th. IT-21 125 (1975) 

A. Chefles, Phys. Rev. A 64, 062305 (2001) 

G. M. D'Ariano, P. Perinotti, and M. F. Sacchi, Euro- 

phys. Lett. 65 165 (2004) 

G. Chiribella, G. M. D'Ariano, P. Perinotti, and M. F. 
Sacchi, Phys. Rev. Lett, in press (quant-ph/0405095 1 
L. Mandel and E. Wolf, Optical coherence and quantum 
optics, Cambridge University Press (Cambridge, New 
York, Melbourne, 1995) 

P. L. Kelley, and W. H. Kleiner, Phys. Rev. 136 A316 
(1964) 

L. Mandel, Proc. Phys. Soc. 72 1037 (1958); ibid. 74 233 
(1959) 

G. M. D'Ariano and P. Lo Presti, Phys. Rev. Lett. 86, 
4195 (2001); 91, 47902 (2003). 

D. N. Klyshko, Sov. J. Quantum Electron. 10, 1112 
(1980) 

A. Migdall, R. Datla, A. Sergienko, J. S. Orszak and Y. 

H. Shih, Appl. Opt. 37, 3455 (1998). 

G. M. D'Ariano, Scuola "E. Fermi" on Experimental 
Quantum Computation and Information, F. De Martini 
and C. Monroe ed. (lOS Press, Amsterdam 2002) pag. 
385. 

K. Banaszek, G. M. D'Ariano, M. G. A. Paris, and M. 

F. Sacchi, Phys. Rev. A 61, R10304 (2000). 
Quantum State Estimation, Lecture Notes in Physics, 
649 M. Paris and J. Rehacek Eds., (Springer, Berlin 
2004) 

This is true in general on theoretical- statistics grounds, 
since when the optimal estimator (i.e. the one achieving 
the Cramer-Rao bound) exists, then it coincides with the 
maximum likelihood estimator. 

G. M. D'Ariano, M. Vasilyev, and P. Kumar Phys. Rev. 
A 58, 636 (1998) 

A. I. Lvovsky and S. A. Babichev, Phys. Rev. A 66, 
011801R (2002) 

J. Wenger, R. Tualle-Brouri, P. Grangier, Opt. Lett. 29, 
1267 (2004) 

A. Zavatta, S. Viciani, M. Bellini eprint 
|quant-p h/0406090 



